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1. Introduction 

The spontaneous chiral symmetry breaking plays a central role in the low-energy Quantum 
Chromodynamics (QCD). It is understood that this phenomena is the source of the hadron masses 
of order Aqcd, the QCD scale. An important exception is the pion, which is nearly massless, 
as it is a pseudo-Nambu-Goldstone boson. The pion dynamics is well described by an effective 
theory, known as chiral perturbation theory (ChPT) ||T]], which is constructed based on the pattern 
of spontaneous symmetry breaking. 

Parameters in ChPT are not known a priori. In phenomenological analysis, they are determined 
with experimental data as inputs, but it is more desirable if they can be calculated starting from the 
first-principles of QCD. This sets a challenge for lattice QCD. At the leading order, there are two 
parameters: the chiral condensate £ and pion decay constant F in the chiral limit. Calculation of 
these parameters has long been one of the main issues in lattice QCD. In particular, the calculation 
of £ has been notoriously difficult, as it survives only in the thermodynamical limit, i.e. the limit 
of massless sea quarks after taking infinite volume limit. The determination of other parameters in 
ChPT, such as the low energy constant at the next-to-leading (NLO) order can be done only after 
the leading order parameters are determined precisely. 

Lattice QCD has become the most powerful tool for non-perturbative calculation of strong 
interaction of hadrons, with the help of the rapid speed-up of computers. In fact, lattice QCD 
has even played a leading role in the development of high-end computers. Still, the mechanism 
of chiral symmetry breaking remained not entirely clear until recently, since the chiral symmetry 
itself was violated in the simulations with the conventional lattice fermion formulations. It is 
theoretically known that the use of Neuberger's overlap-Dirac operator [§] is a solution to this 
problem as it realizes exact chiral symmetry at finite lattice spacings [^,0. Because of its numerical 
cost, however, it was only recently that the large-scale simulation of dynamical overlap fermions 
became feasible. 

The numerical cost of the overlap-Dirac operator is high, compared to other non-chiral or 
non-flavor-symmetric lattice fermions, as it involves an approximation of the sign function of the 
hermitian Wilson-Dirac operator. The cost increases even more when the Atiyah-Singer index of 
the Dirac operator, which corresponds to the topological charge of the background gauge field, 
changes its value by ±1. This is because the molecular dynamics steps have to involve an extra 
procedure in order to catch a sudden jump of the fermion determinant on the topology bound- 
ary. This additional procedure, known as the reflection/refraction, needs numerical cost potentially 
proportional to the lattice volume squared. 

Recently, the JLQCD and TWQCD collaborations have performed large-scale simulations of 
2- and 2-1-1 -flavor QCD employing the overlap fermions for sea quarks We avoid the extra 
numerical cost due to the change of topology by a modification of the lattice action to suppress 
the topology tunneling as proposed in |^ ||, ^. Our lattice simulations are confined in a fixed 
topological sector, so that an expectation value of any operator could be deviated from the value in 
the true QCD vacuum. This effect can be understood as a finite volume effect and estimated in a 
theoretically clean manner as discussed below. It is worth noticed that the simulation parameters 



contain those in the e-regime on a L ~ 2 fm lattice, as well as in the conventional p-regime [10, 11 



12]. This enables us to study the chiral dynamics in an entirely different set-up and to determine 
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the low-energy constants at the point very close to the chiral limit. 

With exact chiral symmetry, the study of spontaneous chiral symmetry breaking is theoretically 
clean, but it still requires a good control of the systematic effects due to the finite volume [13]. For 
such infra-red effects, the lightest particle, which is the pion, gives a dominant contribution. It 
should therefore be possible to use analytic calculations within ChPT in order to predict the finite 
volume corrections for a quantity of interest. Then, the lattice results can be directly fitted with 
these finite-volume formulae of ChPT to determine the relevant low-energy constants. The effect 
of fixed topology can also be understood as one of such infra-red effects since the global topological 
charge should not affect the physics at a local sub-volume when the entire volume V is large enough 
[14, 1^]. In a calculation of the topological susceptibility [16, 17, 18, [l^ through topological 
charge density correlator, we can actually see that local topological excitations are active even when 
the global topological charge is kept fixed. Its result is consistent with an expectation of ChPT, 
which implies that the ChPT-based analysis is valid for the effects due to the fixed topological 
charge[|2|,^. 

There have been a number of analytical works that aimed at controlling the infrared effects 
occurring in the lattice simulations. A well-known example is the finite volume correction due the 
pions wrapping around the lattice [^, 23]. Extended works are necessary when the system enters 
the so-called £-regime [24, 25] by reducing the sea quark mass to the vicinity of the chiral limit. 
In this regime, the vacuum fluctuation of the pion field plays a special role and a non-perturbative 
approach is needed in ChPT. Namely, the zero-momentum pion mode has to be integrated over 
the group manifold of the chiral symmetry in contrast to the case of the conventional /j-regime 
where a certain vacuum is (randomly) chosen by the spontaneous symmetry breaking. Recently, 
the partition functions with fully non-degenerate flavors [26] were calculated, so that even the 
(partially) quenched analysis [ p7| ] of the meson correlators is possible. To study more realistic 
set-up, i.e. including the strange quark in the jc-regime, several hybrid method to treat both the e- 
and /7-regimes have been proposed [28, |^]. The effect of fixed topology is worked out in [21]. We 
also note that the effects of explicit violation of chiral symmetry due to the Wilson term are also 
discussed [30, 31 1, which is needed to study the Wilson fermion simulations near the chiral limit 
[g2|,|33|]. 

In this talk, the dynamical overlap fermion simulation by the JLQCD and TWQCD collabora- 
tions is reviewed in Section ^. In Section ||, we discuss the finite size scaling as well as the global 
topological effects within ChPT. As an example, our recent result for chiral condensate [34, 35, 36] 
is presented in Section ^ Summary and conclusion are given in Section ^. 



2. Dynamical overlap fermion at fixed topology 

We employ the overlap-Dirac operator [^ 

mo + ^j + [mo - -j Y5Sgn[Hw{-mo)], (2.1) 

for the quark action. Here m denotes the quark mass and Hw = YsDwi—mo) is the Hermitian 
Wilson-Dirac operator with a large negative mass —mo. We take mo = 1.6 throughout our simu- 
lations. (Here and in the following the parameters are given in the lattice unit.) In the chiral limit 



3 



Exploring chiral dynamics with overlap fermions 



Hidenori Fukaya 



m — )• 0, the overlap-Dirac operator (2.1) satisfies tlie Ginsparg- Wilson relation 



£»(0)75 + 75£»(0) = — D(0)75£>(0). (2.2) 



With this relation, the fermion action constructed from (2.1) has exact chiral symmetry under a 
modified chiral transformation [Q]. Moreover, it is known that the overlap-Dirac operator has an 
index which corresponds to the topological charge in the continuum limit p7[]. 

In the numerical implementation of the overlap-Dirac operator (p]), the profile of near-zero 
modes of the kernel operator Hw{—mo) largely affects the numerical cost of the overlap fermion 
(The presence of such near-zero modes is also a problem for the locality property of the overlap 



operator [|38|].). For the approximation of the sign function in ([2.1|), the number of operations 



of the Wilson-Dirac operator needed to keep a certain precision monotonically increases as the 
condition number A^^'^/A^'" grows, where x^"^^""" denotes the maximum/minimum eigenvalue 
of the operator \Hw{—mo)\. Moreover, since the overlap-Dirac operator is not uniquely determined 
when Hw{—mQ) has a zero eigenvalue, the overlap fermion determinant has a discontinuity. This 
discontinuity of the determinant prevents smooth evolution of the molecular dynamics steps and 
requires a special treatment, known as the reflection/refraction procedure It needs an extra 
numerical cost, which is potentially proportional to the lattice volume squared. 

At currently available lattice spacings with conventional gauge actions, the spectral density 



Pwi^w) of the operator H\Y(—mQ) is non-zero at zero eigenvalue Aiy = [p9[]. Note that the 



appearance of pwi^w = 0) is, however, a lattice artifact due to the so-called dislocations: local 



lumps of gauge configurations []40|], which disappears in the continuum limit. 

To avoid the problem of the large extra numerical cost and of the potentially ill-defined overlap 
operator, we introduce additional Wilson fermions and twisted-mass bosonic spinors to generate a 
weight 

det[Hwi-mof] 

in the functional integrals [^, ^, Both of fermions and ghosts are unphysical as their masses 
are of order of the lattice cutoff, and thus do not affect low-energy physics. The numerator sup- 
presses the appearance of near-zero modes, while the denominator cancels unwanted effects from 
higher modes. The "twisted-mass" parameter /i controls the value below which the eigenmodes 
are suppressed. In our numerical studies, we set /i = 0.2. 

As Fig. H shows, the near-zero modes of H]v{—m()) are actually washed out when ^ is non-zero 
in quenched QCD simulations. This leads to a large reduction of the numerical cost to approximate 



the sign function in ( [2.1| ) [g]. We also find that the molecular dynamics evolution is smooth in the 
hybrid Monte Carlo updates and we can turn off the reflection/refraction procedure. 

The presence of zero-mode of Hw{—nn)) is related to a topology change: the Atiyah-Singer 
index or the topological chai^ge of gauge fields changes its value when an eigenvalue Hw{—mo) 
crosses zero. The condition Hw{—mo) = 0, thus, forms a topology boundary on the gauge con- 
figuration space. With the lattice action including ( |2.3[ ), therefore, the topological charge never 
changes during the molecular dynamics steps of the Hybrid Monte Carlo (HMC) simulations. In 
this work, the simulations are mainly performed in the trivial topological sector 2 = 0. In order to 
check the topological charge dependence, we also carry out independent simulations at Q = +1, 
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Figure 1: Histogram of the spectral density of Hw{—m(j). Data for three values of ji (ji = 0.0, 0.2, and 0.4) 
are shown in the plot. Note that fi—Q corresponds to the case where the extra fermion determinant is turned 
off. 



—2 and —4 at some parameter choices. The configuration space of a given fixed topology is simply 
connected in the continuum limit, hence it is natural to assume that the ergodicity of the Monte 
Carlo simulation is satisfied within in a given topological sector. 

In the Monte Carlo simulations, we choose 5-6 different points of the up and down quark mass 
Mud in a range 0.002 < m^d < 0. 100. For the Nf = 2 + I runs, two values of the strange quark 
mass: = 0.080 and 0.100 are taken. For the gauge part, we use the Iwasaki gauge action pT] ] at 
j8 = 2.3 (except for the case of m^d = 0.002 in the Nf = 2 run where jS = 2.35 is chosen). The lattice 
volumes are V = L^T = 16^ x 32 {Nf = 2) and V = L^r = 16^ x 48 (A^/ = 2 + 1). For the latter, 
we also carry out a run on a V = L^T = 24^ x 48 lattice at m„j = 0.025 and = 0.080, in order 
to check the finite volume effect. The lattice scales a^^ = 1.667 GeV {Nf = 2) and a^^ = 1.833 
GeV {Nf = 2+1) are determined from the heavy quark potential, using ro = 0.49 fm as an input 
[p2|]. The lattice size is then estimated as L ~ 1.9 fm for Nf = 2, and L ~ 1.7 fm for A/y^ = 2 + 1 
runs. Note that for the Ughtest quark mass m„d = 0.002 ~ 3 MeV, the system of pions is inside the 
e-regime. 

Since our gauge configurations are generated in a fixed topological sector, expectation value 
of any operator could be different from those in the QCD vacuum. Also, our lattice size is ~ 2 fm 
and considerable finite volume effects, especially in the £-regime, are expected. As our lattice size 
is, however, still kept larger than the inverse of QCD scale, i.e. h.QCDL » 1, both effects can be 
considered as a part of infra-red physics for which pions are most responsible. We, therefore, expect 
that chiral perturbation theory (ChPT) can correct these systematic effects. In the next section, we 
discuss how to evaluate physical observables in a fixed topological sector within ChPT at finite 
V . Non-perturbative treatment of the zero momentum mode, as well as the Fourier transform with 
respect to the vacuum angle Q, play a key role. Using their analytic formulae, we can convert the 
lattice QCD results on a finite lattice to the values in the true QCD vacuum in the infinite space-time 
volume. 
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3. Finite V and fixed Q effects within ChPT 

In this section, we first discuss how to evaluate the effect of fixing topology. A general ar- 
gument leads to a consequence that the dependence on the global topological charge only appears 
with a suppression factor l/V . Namely, it is a part of the finite volume effects. Recent studies 
of the finite volume scaling within ChPT are then reviewed. Once we assume that the heavier 
hadrons, such as rho mesons, baryons etc, are all decoupled from the theory at the scale of 1 /V^/'^, 
only pions describe the difference of the finite volume system from the infinite volume one. We 
discuss, in particular, a non-perturbative approach to integrate over the chiral field's vacuum, which 
is necessary in the £-regime. 

3.1 Topology as an infra-red physics 

Let us start our discussion with an intuitively noticeable difference between the trivial topo- 
logical sector {Q = 0) and the first non-trivial one {Q = 1). In the weak coupling limit ^ ^ 1, it is 
well-known that a self-dual solution, the so-called one-instanton solution, dominates the configura- 
tion space of the Q=\ sector and its relative weight is given by ~ exp(— 87r^/g^). For larger value 
of Q, the weight is expected to be ~ exp(— 8;r^|2|/g^). As the coupling constant becomes strong, 
g ~ 1 , more complicated configurations with many pairs of instantons and anti-instantons are more 
favored, since the entropy gives more impact on the free energy than the action density. Suppose 
that the number of such pairs generated in a typical configuration is Qave- The trivial sector 2 = 
then has Qave instantons and Qave anti-instantons while in the 2=1 sector {Qave + 1) instantons 
and Qave instantons are there. As Qave grows, the difference between the global topological charge, 
2 = and 2 = 1 would become less important. 

If the theory has a mass gap Kgap (it is natural to assume Ag^p = ^qcd for the pure gauge 
theory while Kgap is the pion mass ruj^ for QCD), the typical size of an instanton or anti-instanton 
should be given by \/Kgap and their density is estimated as ~ A^ap- The value of Qave discussed 
above is then estimated by ~ ^^gap^ ^^id one can easily see how the difference between 2 = and 
2 = 1 (or higher) disappears as ~ l/V when V is sent to infinity or equivalently Q 

ave — ^ °°- 

The 

effect of the global topological charge thus should be understood as a finite volume effect. 

Brower et al. [14] and Aoki et al. [15] gave a more theoretical and solid formulation for the 
effect of the global topological charge. The partition function of the theory at a fixed topological 
charge Q is obtained from those at the Q vacua by a Fourier transformation 

Zq = j de e'^Qz{e) = j dd e'^Q&xp{-f{e)V), (3.1) 

where f{B) denotes a free-energy density of the 9 vacuum. When the vacuum angle Q is small, 
f{d) can be expanded in as [Q 



= f 0' + C40^ + C60'' + ---, (3.2) 

where a constant term is omitted. Here, Xt corresponds to the topological susceptibihty. Assuming 
that all the constants, Xt, C4, C(, etc. are of the order of ~ (A^ap)^ and the volume is large enough to 
satisfy LAgap ^ 1, the above 6 integral can be evaluated by a saddle-point expansion as 



1 ( 



(3.3) 
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Figure 2: The eta-prime correlator (circles) atm — 0.002 and Q — obtained in the two-flavor QCD simu- 
lation. A negative constant contribution is seen. The triangles are its connected and disconnected diagram 
parts. 



which clearly shows that the global topological charge dependence disappears in the limit V — oo. 
It is also important to notice that the distribution of the global topological charge converges to the 
Gaussian distribution as the volume increases, which agrees well with the intuitive picture above 
that only the entropy given by the distribution of instantons and anti-instantons becomes important 
in the thermodynamical limit. 

Under those minimal assumptions on the vacuum free energy, one can prove that Xt^ ca ^^c. 
can be extracted from lattice QCD simulations at a fixed topological charge For instance, Xt 
appears as a constant mode in the two-point correlator in the flavor singlet channel 

k^^x)ri{y))Q = -| + ^(1/1/2) + ^(^-"vl->i), (3.4) 

for a large separation jj. The excitation in this channel corresponds to the ^' meson whose non- 
zero mass is given by m^t. The constant correlation has a negative sign when the global topological 
charge Q is zero, because at long distances there is more chance to find oppositely charged local 
topological excitations when the sum is constrained to zero. 

In the numerical simulations 1 16, |l7|, [l^, [l^ the presence of this constant mode is confirmed 



as Fig. ^ shows. Moreover, the extracted values of x^ via above formula are found to agree with the 



ChPT prediction [25] 



It = , ' (3.5) 

as seen in Fig. |3[ The value of £ extracted from this analysis is consistent with a nominal value 
£ ~ (250 MeV)^. Chiral fit including the next-to-leading chiral corrections [20, ^ is underway. 

There are two remarkable conclusions that may be drawn from these lattice data. First, lo- 
cal fluctuation of topology exists even when the global topological charge is fixed in Monte Carlo 
simulations. There was some doubt about the ergodicity of the Monte Carlo simulation with the 
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Figure 3: Xt extracted from the tj' meson correlators. A good agreement with ChPT predictions (solid lines) 
is seen both in the A^/ = 2 and Nf = 2 + \ lattice data. 



topology fixing term, but as far as tlie numerical data imply there is no evidence of the problem. 
Second, the topological charge actually feels the presence of dynamical fermions and the Xt van- 
ishes in the chiral limit as expected from ChPT. Topology is a part of the infrared physics that can 
be well described by the pion physics. 

3.2 Finite V and fixed Q within ChPT 

The Lagrangian of ChPT is given by 

^ = ^Tv[d^U{xfd^U{x)] - ^Tr[^^e-'^l^fU{x) + U{x)'^ e'^'^f + • • • , (3.6) 

where the chiral field U (x) is an element of SU {Nf) group. Here the pion decay constant and the 
chiral condensate are denoted by F and Z, respectively. The vacuum angle 6 is given as a phase in 
front of the mass matrix ^ = diag(m„, mrf,mv, • • •). 

In the conventional /^-expansion, we treat the exponent of U (x) as the Nambu-Goldstone 
modes (here we denote as i§ (x)), or pions, and expand the chiral field as 

t/(x) =exp(^/^^^j =l + /^^(x)--1^2(x) + .... (3.7) 

With the counting rule 

J^^p\ d^^p, l/L,l/r~p, ^{x)^p, (3.8) 

physical amplitudes are systematically expanded in terms of p^. 

In the p-regime, the finite volume effect appears in the pion propagator since the momentum 



space is discretized [22|. Pion correlator reads 



Jp(x-y) 

(§«(x)^^3')) = 5«.£^— 3, (3.9) 
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where a{b) denotes the a{b)-th generator of SU (Nf) and the summation is taken over the 4-momentum 
p = 27i{nt/T,n^/L,ny/L,n^/L), with integer n^'s. As a consequence, all the correlators become 
periodic. Even at a contact point x = y, there exists a finite volume correction 

= l4^^i(^l«l)' (3.11) 

which is understood as an effect of pion wrapping around the lattice. Here (x) is the modified 
Bessel function and the summation is taken over the 4-vector = n^L^ with L, = L for / = 1,2,3 
and L4 = T, except for = (0,0,0,0). Note that the subtraction of the ultraviolet divergence is 
done at a scale Hsub, which can be made in exactly the same way as in the infinite volume. In 
a similar perturbative manner, the effect of global topology is recently calculated to the next-to- 
leading order ||2l|]. 

In the £-regime, the above p-expansion ( ^^ ) fails because the zero-momentum mode contri- 
bution induces an unphysical infrared divergence, which has to be circumvented by exactly treating 
the vacuum fluctuation of the chiral field. Namely, using a parameterization 

U{x) = [/oexp I y ' , (3.12) 

where Uq G SU (Nf) and satisfies 

J d^x^'{x) = 0, (3.13) 

one can explicitly factorize the zero momentum part as Uq. Since Uq has no dependence on x, the 
group integral can be non-perturbatively performed as in the calculation of random matrix models. 
The non-zero momentum modes i^"s are perturbatively treated as an expansion in according to 
the counting rule 



^4 



e\ d^r^e, l/L,l/r~e, (^'(;c)~e. (3.14) 



This £-expansion [gj, Q is useful when the quark mass is so small that the pion correlation length 
exceeds the spatial extent, mj^L <^ 1 . 

The zero momentum component Uq can be explicitly integrated out and written in terms of 
analytic functions. This fact opens an interesting theoretical opportunities. In particular, at the 
leading order of the £-expansion, the system is proven to be equivalent to the Random Matrix 
Theory [44, In the context of the QCD study, this provides a new method to determine the 

chiral condensate by matching the low-lying eigenvalues of the Dirac operator. At an early stage, a 
simple setup with all degenerate quark masses were studied. As lattice QCD is developed to reach 
the simulations near the chiral limit, calculations in a more realistic setup has become relevant, and 



partially quenched calculations of various quantities have been carried out [ p7| , With strange 
quark mass kept at its physical value, the finite volume system is not purely in the £-regime even 
when the up and down quark masses are sent close to the chiral limit, because the kaon and 77 are 
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heavy and do not satisfy mK,nL <^ 1. For this mixed-regime, a hybrid method is proposed [28] and 



even extended to the case of heavy-light mesons [148|]. More recently, a theoretical framework in 
which the e- and jc-regimes are treated in a unified manner is proposed [ |2^ ] of which details are 
described in the next section. 

As a final remark of this section, we note the role of topological charge in the e-regime. 
Intuitively, the global topological charge become relevant to the dynamics of the system when the 
volume is small. This can be explicitly studied within ChPT. The 6 integral can be absorbed in the 
zero-mode integrals 

f^f dUoe'^Q=f dUoidetUof, (3.15) 

J 271 JsU{Nf) Ju(Nf) 

where Uq is integrated over \J{Nf) manifold. The effect of the topological charge enters through 
a factor {detUo)^. For instance, the spectral density of low-lying Dirac eigenmodes is largely 
affected by Q, of which dependence can be used to test the validity of ChPT, in addition to the 
quark mass dependence. 

4. Determination of the chiral condensate 

4.1 Analytic results beyond the leading order 

The chiral condensate is related to the Dirac eigenvalue density p(A) at A = in the thermo- 
dynamical limit [ p9| ] as p(0) = 11/71. This relation can be easily extended to non-zero eigenvalues 
by an analytical continuation of the valence mass m,, to a pure imaginary value iX : 

p(A) = ^Re(^-,^„>|,n„=,A- (4.1) 

Here, qvqv is the scalar density operator made of the valence quark field. This general formula is 
valid for both p- and £-regimes. 

In the /^-regime, using the partial quenching technique for the imaginary valence quark mass, 
Osborn et al. [50| (see also [^T]]) found that the Dirac spectrum contains a logarithmic dependence 



on A. This calculation is done in the infinite volume limit with degenerate quark masses. 

For small eigenvalues, the effect of finite volume becomes important. The ChPT calculation 
is simplified if one consider the e-expansion and taking its leading order contribution. The integral 
over the zero momentum pion mode can be done analytically, and the spectral function has been 
obtained as a function of Nf, sea quark masses, and topological charge Q ||5^, 53, |^]. Except for 



the exact zero-modes associated with Q, there is a finite gap from zero (of order l/'LV , which is 
called the microscopic region) in the Dirac operator spectrum. These analytic ChPT results can 
be used to extract Z by comparing with the lattice data in the e-regime 13 1. But, since the 



formulae are obtained at the leading order, the value of £ thus obtained is a subject of the NLO 
corrections of the £-expansion. Furthermore, it requires that the system is in the e-regime, which is 
numerically demanding. For common lattice QCD configurations produced in a /j-regime set-up, 
these analytical results cannot be applied. 



Here we introduce a new method of the chiral expansion [29]. It is based on the ;7-expansion. 



but includes the pion zero-mode integral explicitly so that a transition to the e-regime is smooth. 
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In this scheme, one may predict the eigenvalue spectrum in the microscopic region for the system 
in the p-regime. With the so-called replica trick, the calculation is extended to the case of non- 
degenerate quarks of arbitrary number of flavors. 
At a fixed topological charge Q, we obtain [ p9| ] 

Pe(A) = £effPe(Ai:effV,{m,,„reftV}) + p^(A,{m,,4), (4.2) 

where A denotes the Dirac eigenvalue, {m^eaZeff^} = {w^i^eff^;'«2^eff^5 • ■ "} is a set of the sea 
quark masses normalized by an effective chiral condensate Zgff (the definition is given below) and 
V. The first term on the right hand side of (1.2) contains the one in the leading-order £-expansion 
Pq(C) {Mseo} = {Mi5M2 ■■■})> which is rescaled so that the physical scale £ is factored out. This is 
a known function given by determinants of the Bessel functions |54]: 



Pe(?,{M.«J) = C2- 



det,: 



where an Nf x Nf matrix £/ and an {Nf + 2) x (Nf + 2) matrix are defined by 



(4.3) 



'2j 



(//1,2). 



(4.4) 



(4.5) 



Here, the overall sign is C2 = +1 for the Nf = 2 and 3 cases. 

The second term in ( |4.2| ) is a logarithmic NLO correction as always seen in the conventional 
p-expansion. Defining M?- = (m,- + mj)L/F^, the function is given by 



P^(A,{m,,4) 



Nf 

Re[£(A(M- 

/ 



-A{Mjf/2))-{G{Ml 



-G(0))] 



(4.6) 



where 



2 
3 



A(m2) + (M- 

2(M2,-M2^2 



+ 



{M^-Mlj){M-^-Ml) 



A(M2) + 

dM2A{M^ 



9(A/2_^2)2 



16712 



M 



A(M2 



(A^/ = 2), 



(^/ = 2 + l), 



(4.7) 
(4.8) 



Here M^^ = 2m,,L/F^- = 2mj£/F2, m2 = 2m,L/F^- and = (M^^ + 2M2)/3. The function 
gi (M^) = gi (M^) — 1 jM^V denotes the well-known finite volume correction from non-zero modes 
[ p2| ] (see also ( |3.10| )). The scale ^^^^(=770 MeV in this work) is a subtraction scale. Note that 
p''(A, {mvea} is insensitive to the topological charge. 

The effective condensate in (4^) is also expressed in terms of ~A(M^) and G(M^^ as 



-eff 



1 



1 - ( £A(m2^/2) - (7(0) - \^U,Y^M 



ff 



(4.9) 
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Figure 4: The first term EgffjQg (solid-thin curve), the second term pP (dashed) and the total contribution 
(solid-thick) of the spectral density ( 4.2 ) are shown. The curves are multiplied by n. We use E — [240MeV]^, 
F = 94 MeV, LI = -0.0001, L = r/3 = 1 .9 fm, = 20 MeV and m, = 120 MeV as inputs. 



This depends on the sea quark masses, volume V and Lg of which value is renormalized (at }Xsub = 
110 MeV in this work). 

For an illustration, we draw curves given by the formula (O) in Fig. Q The contributions 



from the first term ZeffPg (sohd-thin curve), the second term pP (dashed) and the total contribution 
Pq(A) (solid-thick) are shown separately. We use typical parameters £ = [240MeV]^, F = 94 MeV, 
LI = -0.0001, L = r/2 = 1.9 fm, niud = 20 MeV and m,, = 120 MeV as inputs. One can see that 
the second term gives a negative contribution and shows a significant curvature in the lower end of 
the spectrum. This is the effect of the chiral logarithm. For this quark mass, the formula starts to 
deviate from the leading order expression in the £-expansion already at A ~ 5 MeV. 

4.2 A numerical analysis 

Our simulation details and parameters have been already presented in Section. |2[ For the 
study of the Dirac spectrum, 80 lowest pairs of eigenvalues of the overlap-Dirac operator D{0) are 
calculated at every 5-10 trajectories. We employ the implicitly restarted Lanczos algorithm for 
the chirally projected operator /'+D(0)P+, where P^ = {1 + Ys)/!. From its eigenvalue ReA°^', the 
pair of eigenvalues X"^' (and its complex conjugate) of D{0) is extracted through the relation |1 — 
A'^' /mop = 1, that forms a circle on a complex plane. For the comparison with the effective theory, 
the lattice eigenvalue is projected onto the imaginary axis as A = 1mA "^'/(l — ReA''^'/(2mo)). 
Note that the real part of A"^' is negligible (within 1%) for the low-lying modes. 



When we match the lattice data for the spectral density with the analytic calculation (4.2), 



two parameters are to be determined at each set of the quark masses: Zeff and F. In the second 



NLO term of (1.2), the difference between Zeff and £ is a higher order effect. We therefore take 
two reference values of A to give inputs to determine Zeff and F. The reference points are chosen 
such that they have maximum sensitivity to the parameters in the convergence range of the chiral 
expansion: A = 0.004 (~ 7 MeV) and 0.017 (~ 30 MeV) except for the case with niud = 0.002 
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and 2 = 1, for which we choose A = 0.01 and 0.02 (because of its weaker sensitivity to the NLO 
effects). At these two reference points, we compare the mode number below a given value of A 



[g3p, with an integrated formula of ChPT 



Nq{X) = V / dX'pQ{X'), (4.10) 

JO 



and determine Zgff and F . As Giusti and Liischer [ p3[ ] studied, it is also useful to define a quantity 



^ (4.11) 

to see the NLO effects, or the chiral logarithmic effects to £. We test the both of A^/ = 2 + 1 and 
Nf = 2 ChPT formulae. For the latter case, the strange quark is assumed to be decoupled from the 
theory. 

Figures || and ^ show the lattice data for the spectral density (upper panel), its integral (middle) 
and 1!q"^''{X) defined by ( [4.1 1 ) at two different sea quark masses: one in the /j-regime (m = 0.015, 



Fig. Ill) and the other in the £-regime {m = 0.002, Fig. The analytic formula is also plotted with 
two parameters fixed at two reference points of the mode number. The leading-order contribution 
is given by dotted curves while the full result is shown by solid curves. 

In the ;7-regime result (Fig. the effect of the NLO term in the /^-expansion is clearly seen 
as a deviation from the leading-order density TeffPg (dotted curve) in the histogram. The deviation 
starting already around A ~ 0.005 is also clear in the mode number Nq{X) and Y!q'^''{X). On the 
other hand, the NLO formula (solid curve) describes the lattice data very nicely up to A ~ nis/2. 

The convergence of the chiral expansion is better for the £-regime data (Fig. ^, but the differ- 
ence between LO and NLO still exists. We also observe that there is a wider gap near A = 0, which 
is expected because the value of the sea quark mass m = 0.002 is similar to the lowest eigenvalue, 
so that the suppression due to the fermionic determinant Y\{X^ + m^) works strongly. 

One of the significant consequences of the ChPT formula ( |4.2[ ) is that the spectral function for 
different topological charge Q and volume V should be described by the same set of the parameters, 
i.e. Zeff and F. This provides a highly non-trivial cross-check of the formula. For this purpose we 
produced data at non-zero topological charge 2 = 1. The results are shown in Fig. ^ Here the 
curves of the NLO ChPT is drawn with inputs from the Q = data and there is no further free 
parameter to adjust. The good agreement below A ~ 0.03 gives further confidence on the analysis. 

A similar check can be done with the lattice data obtained from a larger volume lattice 24^ x 
48, for which the data are shown in Fig. |8[ The comparison is a bit more tricky for different 
volumes, because the definition of Zeff ( |4-9| ) depends on V . Namely the function A(M^) contains 
gi(M^), which represents the finite volume effect. It is possible to convert the value of Zeff for 
different volumes. If we convert the result at L = 24, Zgff = 0.00306(7) to the one on a L = 16 
lattice, it becomes 0.00341(18), which may be compared with the independent calculation at L = 16 
at the same sea quark mass m = 0.025, which is 0.0333(18). Therefore, the finite volume scaling 
is confirmed at least on two different volumes, whose difference is a factor of 3. 

The curves in Figures Q-^ are drawn using the Nf = 2 + \ ChPT results, but we found the 
difference from Nf = 2 ChPT formula is hardly visible in the scale of this plot, which confirms 
decoupling of the strange quark from the low energy theory. 
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a 0.0015 



0.06 




0.06 



Figure 5: The spectral density %Pq{X) (top), the mode number A^q(A) (center) and Eg"''*(A) (bottom) 
of the Dirac operator at niud = 0.015, nis = 0.080 and 2 = 0. The lattice result (histogram (top) or solid 
symbols (center and bottom)) is compared with the ChPT formula drawn by solid curves. For comparison, 
the prediction of the leading e-expansion (dashed curves) is also shown. 
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Q. 



0.005 



0.004 



0.003 



0.002 



0.001 




0.01 0.02 0.03 0.04 0.05 0.06 




a 0.0015 



0.06 



0.01 0.02 0.03 0.04 0.05 0.06 



Figure 6: Same as Fig. ^ but at = 0.002. The NLO correction is smaller in the e-regime. 
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Q. 



0.005 



0.004 



0.003 



0.002 



0.001 



0.0015 



lattice 
ChPT LO 
ChPTNLO 




yn^=0.Q\5, e=0 
ni„^=0.015, e=l 
ChPT NLO 



0.02 0.03 0.04 0.05 0.06 

A, 




0.01 0.02 0.03 0.04 0.05 0.06 



Figure 7: The spectral density at m = 0.015 and 2=1 (top) and comparison of Y!^'^'^{X) at Q = and 1 
(bottom). In the ChPT curves, the same values of E^^j and F are used as inputs. 



From these analysis the values of Zeff and F ai^e extracted for each sea quark mass. Note that 
Zeff is extracted at the NLO accuracy, while the value of F, which first appears at the NLO term, 
might have larger systematic corrections from NNLO contributions. We find that the results for 
Zeff are stable under change of two reference points in a range A < 0.03. As noted above, there is 
little difference between Nf = 2 and Nf = 2 + \ formulae; Zeff and F are almost equal well within 
the statistical error. The difference between = 0.080 and rris = 0.100 is even weaker. In the 
following analysis, we concentrate on the data at rus = 0.080. 



4.3 Chiral extrapolation of Zeff 



We next consider the sea quark mass dependence of Zeff- As shown in (19), Zeff is a function 
of r, F, Lg, which can be determined from the lattice data. The chiral condensate £ thus obtained 
should have the NLO accuracy. In the fitting of the lattice data, we attempt (A) 3-parameter (£, F, 
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Q. 



0.005 



0.004 



0.003 



0.002 



0.001 



0.0035 



lattice 
ChPT LO 
ChPTNLO 




0.01 0.02 0.03 0.04 0.05 0.06 

A, 




0.06 



Figure 8: The spectral density at m = 0.025 and L = 24 (top) and comparison of TJ^'^''{K) at L = 24 and 
16 (bottom). In the ChPT curves, the same values of T^^ff (but the volume dependence is corrected within 
ChPT) and F ai^e used as inputs. 



Le) fit without any inputs and (B) 2-parameter (I, Le) fit with F = 0.0410 (for A'/ = 3 ChPT) or 
with F = 0.0406 (for Nf = 2 ChPT). These values of F correspond to the chiral limit of F extracted 
from the analysis of the spectral function. 

The fitting is shown in Fig. ^ for the case (A) with the Nf = 3 ChPT formula. We use the 
lightest 4, 5, and 6 data points. All the curves are consistent with the lattice data used in the fit and 
in fact the per degrees of freedom is reasonable (between 0.6 and 1.5). A remarkable fact is that 
the chiral limit (shown by a square) is not sensitive to the number of data points used. The chiral 
limit is very stable because of the presence of the £-regime data point. Similar curves are obtained 
for the case with Nf = 2 and for the case (B). With the 2-parameter fit (the case(B)) the heaviest 
data point cannot be well described, i.e. ;^^/d.o.f. is about 2.5. 

From these curves, one can extract the low energy constants of ChPT. Note in the case of 
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0.006 r- 




Figure 9: Three parameter fit of Eeff . The Nf = 3 ChPT formula is used. 

Nf = 3 ChPT, we have two limits of chiral condensate: r^/=^, where "three" flavor massless limit 
is taken, and EPi^y^, which is a two-flavor chiral limit with strange quark mass fixed at a finite value 
Ms = 0.08. As already mentioned, the strange quark dependence is so small that the difference 
from the value at the physical strange quark mass is negligible. The extracted values of Z^^'^^ are 
stable against the different choice of fitting function and fitting range, while L^f=^ shows strong 
sensitivity to them. It means that the determination of L^f=^ is not feasible with our current data 
set. This is natural because the strange quark mass dependence is not well controlled by the lattice 
data. On the other hand, the determination of L^^'^^ is very stable, thanks to the £-regime data point. 
Our estimate of systematic effects due to the chiral extrapolation is ~ 2 %. 

From the above analysis, we determine the low-energy constants for 2-1-1 -flavor QCD as 

£phys ^ 0.001 86(10) (44) ~ [226(4)(18)MeV]^ (4.12) 
F = 0.0406(05)(41) ~ 74(l)(8)MeV, (4.13) 
L^(770MeV) = -0.0001 1(25) (11), (4.14) 

where the first error is statistical and the second error is systematic, respectively. 

To obtain the final result, we convert the value of Z^^'^^ to the definition in the MS scheme. 



by using the non-perturbative renormalization factor [56|calculated through the RI/MOM scheme 



The result [34], IP^y^ in the limit of m„j = and nig fixed at its physical value, is 

£'^(2 GeV) = [242(04)(+}^) MeV]^ (4.15) 



Let us here discuss possible systematic errors in ( |4.15| ). Since our lattice studies are done 
at only one value of j3, it is difficult to estimate the discretization errors. But it should be partly 
reflected in the mismatch of the observables measured in different ways. We here estimate it from 



a mismatch of the lattice spacing; 0.1003(46) fm from the pion decay constant [57] and 0.1087(15) 



fm from the Q. baryon mass [58]. This 7.4% deviation is added in the systematic error. The 
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systematic error due to finite volume is estimated as ~ 1 .4% using the lattice data at two different 
volumes. 

5. Summary and Conclusion 

In this talk, a study of the spontaneous chiral symmetry breaking performed by the JLQCD and 
TWQCD collaborations has been presented. We discussed that fixing topology is an essential part 
for the dynamical overlap fermion simulations on the lattice. By reducing the numerical cost with 
the topology fixing determinant, we have performed the first large-scale simulations of dynamical 
overlap quarks. The up and down quark masses are reduced close to the physical point. We have 
then discussed that the global topology, as well as the finite volume effect, can be described well 
within the chiral perturbation theory. In fact, we have found a good agreement of our lattice data 
for the Dirac operator spectrum with the ChPT predictions even in the region where its finite size 
effect is large. We extract the chiral condensate in 2+1-flavor QCD. 
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